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Abstract 

A systematic algorithm for building integrating factors of the form fj,(x, y), [i(x, y') or y') 
for second order ODEs is presented. The algorithm can determine the existence and explicit 
form of the integrating factors themselves without solving any differential equations, except 
for a linear ODE in one subcase of the fi(x,y) problem. Examples of ODEs not having point 
symmetries are shown to be solvable using this algorithm. The scheme was implemented in 
Maple, in the framework of the ODEtools package and its ODE-solver. A comparison between 
this implementation and other computer algebra ODE-solvers in tackling non-linear examples 
from Kamke's book is shown. 

1 Introduction 



Although in principle it is always possible to determine whether a given ODE is exact (a total 
derivative), there is no known method which is always successful in making arbitrary ODEs exact. 
For n th order ODEs - as in the case of symmetries - integrating factors (/x) are determined as 
solutions of an n th order linear PDE in n+1 variables, and to solve this determining PDE is a 
major problem in itself. 

Despite the fact that the determining PDE for [i naturally splits into a PDE system, the problem 
IS - clS cl whole - too general, and to solve it a restriction of the problem in the form of a more 
concrete ansatz for \i is required. For example, in a recent work by 0] the authors explore possible 
ansatzes depending on the given ODE, which are useful when this ODE has known symmetries of 
certain type. In another work, Q explores the use of computer algebra to try various ansatzes for 
/i, no matter the ODE input, but successively increasing the order of the derivatives (up to the 
n th —1 order) on which ji depends; the idea is to try to maximize the splitting so as to increase 
the chances of solving the resulting PDE system by first simplifying it using differential Groebner 
basis techniques. 



Bearing this in mind, this paper presents a method, for second order explicit ODEs[], which 
systematically determines the existence and the explicit form of integrating factors when they 
depend on only two variables, that is: when they are of the form /i(x,y), n(x,y') or fi(y,y'). The 
approach works without solving any auxiliary differential equations - except for a linear ODE in one 
subcase of the (i(x, y) problem - and is based on the use of the forms of the ODE families admitting 
such integrating factors. It turns out that with this restriction - /i depends on only two variables 
- the use of differential Groebner basis techniques is not necessary; these integrating factors, when 
they exist, can be given directly by identifying the input ODE as a member of one of various related 
ODE families. 

The exposition is organized as follows. In sec. the standard formulation of the determination 
of integrating factors is briefly reviewed and the method we used for obtaining the aforementioned 
integrating factors fj,(x,y), [i(x, y') or y') is presented. In sec. ||, some aspects of the integrating 
factor and symmetry approaches are discussed, and their complementariness is illustrated with 
two ODE families not having point symmetries. Sec. || contains some statistics concerning the 
new solving method and the second order non-linear ODEs found in Kamke's book, as well as a 
comparison of performances of some popular computer algebra packages in solving a related subset 
of these ODEs. Finally, the conclusions contain some general remarks about the work. 

Aside from this, in the Appendix, a table containing extra information concerning integrating 
factors for some of Kamke's ODEs is presented. 

2 Integrating Factors and ODE patterns 

In this paper we use the term "integrating factor" in connection with the explicit form of an 
n th order ODE 

j/W -$(x,y,y>,...,y( n -V) = (1) 
so that /i(x,y,y', ...,y ( - n ~ 1 ^) is an integrating factor if 

n(yM-$) = -^R(x,y,y',...,y^) (2) 

for some function R. The knowledge of fi is - in principle - enough to determine R by using standard 
formulas (see for instance Murphy's book). To determine fx, one can try to solve for it in the 
exactness condition, obtained applying Euler's operator to the total derivative H = fi (y^ — &j : 

dH d (dH\ d 2 (dH\ d n / dH \ _ 

~dy~~dx~ (jh/J + W7 + "' + dx^ [dy^J ~ ( } 

Eq. (^) is of the form 

A(x,y,y',...,y( 2n - V) + y^ n ^B(x, y, y' , y^) = (4) 

1 We say that a second order ODE is in explicit form when it appears as y" — <&{x,y,y') = 0. Also, we exclude 
from the discussion the case of a linear ODE and an integrating factor of the form /-i(x), already known to be the 
solution to the adjoint ODE. 
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where A is of degree n — 1 in y^ and linear in y^ for n < k < (2n — 3), so that Eq.(Q) can be 
split into a PDE system for \i. In the case of second order ODEs - the subject of this work - Eq.(||) 
is of the form 

A(x,y,y')+y" B(x,y,y') =0 (5) 
and the PDE system is obtained by taking A and B equal to zercQ: 

A = (y'fiy'y - fiy + [lyi x ) $ + (®y' x + Ij'^y'y ~ ®y) fJ, + I/' 2 fJ-yy 

+ (fJ-y^y' + fJ-y'^y + 2 fl x y) y' + fJ,y'^ x + [J, x ® y > + fJ, xx = (6) 
B = y'fJ-y'y + $ /iy'j,' + /J- $yly> + 2 fJLy + 2 + = (7) 

Regarding the solvability of these equations, unless a more concrete ansatz for fj,(x, y, y') is given, 
the problem is in principle as difficult as solving the original ODE. We then studied the solution 
for /i of Eqs.(||) and (0) when /i depends only on two variables, that is: for and 
/j,(y,y'). Concretely, we searched for the existence conditions for such integrating factors, expressed 
as a set of equations in <]?, plus an algebraic expression for \i as a function of <!>, valid when the 
existence conditions hold. Formulating the problem in that manner and taking into account the 
integrability conditions of the system, Eqs.@ and (0) turned out to be solvable for fj,(x,y), but 
appeared to us untractable when /i depends on two variables one of which is y' . 

We then considered a different approach, taking into account from the beginning the form of 
the ODE family admitting a given integrating factor. As shown in the following sections, it turns 
out that, using that piece of information (Eq.(|ll|) below), when /i depends only on two variables 
the existence conditions and the integrating factors themselves can be systematically determined; 
and in the cases n(x,y') and fi(y,y'), this can be done without solving any differential equations. 

Concerning the ODE families admitting given integrating factors, we note that, from Eq.(||) 

,(x,y,y',...,y^) = ^ (8) 
and hence the first integral R is of the form 

R = G{x,y,...J n ~ 2 ))+ J fidy^ (9) 
for some function G. In turn, since R is a first integral, it satisfies 

R x + y , Ry+... + $R y (n-l) =0 (10) 

Inserting Eq.(^) into the above and solving for y^ leads to the general form of an ODE admitting 
a given integrating factor: 

In a recent work by |2[, the authors arrive at Eq.(|) and Eq.(§) departing from the adjoint linearized system 
corresponding to a given ODE; the possible splitting of Eq.(^) into an overdetermined system for fi is also mentioned. 
However, in that work, y" of Eq.(^|) above appears replaced by $(a;, y, y'), and the authors discuss possible alternatives 
to tackle Eqs.(||) and (fjj) instead of Eqs.(|) and (0). 
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(ii) 



2.1 Second order ODEs and integrating factors of the form fi(x,y) 

We consider first the determination of integrating factors of the form /i(x,y), which turns out to 
be straightforward^. The determining equations (|6|) and (|^) for this case are given by: 

y' 2 Hyy + 2 H xy y' + fl ®y< x + (1 ®y> y y' - fJ,$y- fJ.y$ + [ly^ $ <y> + fJ. xx + (J, x ®y> = 

(12) 

/i §ylyl +2fly = 



Although the use of integrability conditions is enough to tackle this problem, the solving of Eqs.(12) 
can be directly simplified if we take into account the ODE family admitting an integrating factor 
[i(x,y). From Eq.(|Tl|), that ODE family takes the form 



y" = a(x, y) y' 2 + b(x, y) y' + c(x, y), (13) 



where 



a(x,y) = -^, b (x,y) = -^±^, c (x,y) = (14) 
fx, fj, fx 

and G(x,y) is an arbitrary function of its arguments. Hence, as a shortcut to solving Eqs.(|i~2|). 
one can take Eq.(O) as an existence condition - $ must be a polynomial of degree two in y' - and 
directly solve Eqs.(|l4|) for fi. The calculations are straightforward; there are two different cases. 

Case A : 2a x — b y ^ 

Defining the two auxiliary quantities 

tf = c y -ac-b x , T = a xx + a x b + <p y (15) 
an integrating factor of the form fi(x, y) exists only when 

T y -a x = 0, T x + ip + bT -T 2 = (16) 
and is then given in solved form, in terms of a, b and c by 

fj,(x,y) = exp (^J ^-T + J adyj dx - J ady^j (17) 

So, in this case, when an integrating factor of this type exists there is only oneQ and it can be 
determined without solving any differential equations. 

Case B : 2a x — b y = 



3 The result for Case A presented in this subsection is also presented as lemma 3.8 in ||. 

4 We recall that if fi is an integrating factor leading to a first integral tp, then the product fj,F(ip) - where F is an 
arbitrary function - is also an integrating factor, which however does not lead to a first integral independent of ip. 
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Redefining ip = c y — a c, an integrating factor of the form fx(x, y) exists only when 

a xx - a x b - = 0, (18) 

and then fi(x,y) is given by 



a dy 

H(x,y) = v(x) e J (19) 
where v(x) is either one of the independent solutions of the second order linear ODE[] 

v" = A(x) v' + B(x) u, (20) 

and 

A(x) = 21-b, B(x) = <p+(l--^\{b-X), 2=^J ady (21) 

So in this case, to transform Eq.([i~9|) into an explicit expression for n we first need to solve a 



second order linear ODE. When the attempt to solve Eq.(20) is successful, using each of its two 
independent solutions as integrating factors leads to the general solution of Eq.([0|), instead of just 
a reduction of order. 



2.2 Second order ODEs and integrating factors of the form fi(x, y') 

When the integrating factor is of the form /i(x,?/), the determining equations (0) and (0) become 

(22) 

<£> jJLylyl + /i Qyly> + 2 fly'^y' + = 

As in the case fi(x,y), the solution we are interested in is an expression for jj,{x,y') in terms 
of as well as existence conditions for such an integrating factor expressed as equations in <3?. 
However, differently than the case fj,(x,y), we didn't find a way to solve the /i(x,y') problem 
just using integrability conditions, neither working by hand nor using the specialized computer 
algebra packages diffalg Q and standard-form Q . We then considered approaching the problem as 
explained in the previous subsection, departing from the form of Eq.(|9|) for fi = fi(x,y'): 

a i i\ F x + Gx + Gy y' . . 

y" = $(x, y, y') = y — (23) 

where G{x, y) and F(x, y') are arbitrary functions of their arguments and 

»(x,y')=F y i (24) 

Now, Eq.(p3[) is not polynomial in either x, y or y', and hence its use to simplify and solve the 
problem is less straightforward than in the case n(x,y). However, in Eq.(|2^), all the dependence 
on y comes from G(x, y) in the numerator, and as it is shown below, this fact is a key to solving 

5 When the given ODE is linear, Eq.(EC]) is just the corresponding adjoint equation. 
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the problem. Considering ODEs for which § y ^ 00, the approach we used can be summarized in 
the following three lemmas whose proofs are developed separately for convenience. 



Lemma 1. For all linear ODEs of the family Eq.(28), an integrating factor of the form n(x,y') 
such that fiyi ^ 0, when it exists, can be determined directly from the coefficient of y in the input 
ODE. 



Lemma 2. For all non-linear ODEs of Eq.(2S), the knowledge of fi(x,y') up to a factor depending 
on x, that is, o/^ r (x,y / ) satisfying 

Hx , y >) = ^yl (25 ) 

is enough to determine H (x) by means of an integral. 



Lemma 3. For all non-linear ODEs members of Eq.(25), it is always possible to determine a 



function J-(x,y') satisfying Eq.(2t) 



Corollary. For all second order ODEs such that & y ^ 0, the determination of fi(x,y') (if the ODE 
is linear we assume fi y ' ^ 0), when it exists, can be performed systematically and without solving 
any differential equations. 

2.2.1 Proof of Lemma 1 

For Eq.(^) to be linear and not missing y, either G x or G y must be linear in y. Both G x and G y 
cannot simultaneously be linear in y since, in such a case, G x /F y i or y'G y /F y i would be non-linear 
in {y, therefore, either G yy = or G xy = 0. 

Case A : G y is linear in y and G xy = 

Hence, G is given by 

G = C 2 y 2 + dy + g(x) (26) 



where g(x) is arbitrary. From Eq. (f23D , in order to have y'G y /F y i linear in {y,y'}, F y > must of 
the form v{x)y' for some function v(x). Also, F x /F y i can have a term linear in y', and a term 
proportional to 1/y' to cancel with the one coming from G x /F y / = g'/F y >, so that 



/ /2 

F y > = uy', F x = -^-g' (27) 

where the coefficient v' /2 in the second equation above arises from the integrability conditions 
between both equations. Eq.(|23|) is then of the form 

y = -7T- y y 28 

2 v v v 



6 ODEs missing y may also have integrating factors of the form /i(x,y'). Such an ODE however can always be 
reduced to first order by a change of variables, so that the determination of a fi(x, y') for it is equivalent to solving a 
first order ODE problem - not the focus of this work. 

7 We are only interested in the case fj, y i = F y i y i 7^ 0. 
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and hence, a linear ODE y" = a{x)y' + b{x)y has an integrating factor fi(x, y') = y'/b when 
b'/b - 2a = 0. A 

Case B : G x is linear in y and G yy = 

In this case, in order to have Eq.(|23|) linear, F y > cannot depend on y', so that the integrating factor 
is of the form fj,(x) and hence the case is of no interest: we end up with the standard search for 
n{x) as the solution to the adjoint of the original linear ODE. A 

2.2.2 Proof of Lemma 2 

It follows from Eqs.(p3|) and (p4|) that, given T satisfying (12 



£ f *{ XiV j) h*A = - G v*( x >y)+ G w(*>yw (29) 

Hence, by taking coefficients of y' in the above, 

W = $ v (x, y, y') F{x, y>) -y'-^ U y {x, y, y') F{x, y')\ = -°j4^A 

V2 = ®L y{x ,y,y>)F^y>)\=- G yy^ (30) 

where the left-hand-sides can be calculated explicitly since they depend only on $ and the given 
T . Similarly, 

d ,,\ F y , x (x,y') + G y (x,y) 
if4 = —J : (x,y) = — - 31 

Q y n(x) 

Now, since in this case the ODE family Eq.(p3[) is nonlinear by hypothesis, either ip2 or 994 is 
different from zero, so that at least one of the pairs {/pi, (^2} or {^3, (^4} can be used to determine 
H(x) as the solution of a first order linear ODE. For example, if ip% 7^ 0, 

— (ipi(x,y) = — (Mx,y) M(x)) (32) 

from where 

1 f^pi _ <teA dx 

., <p 2 \dy dx J 
= e (33) 

If (^2 = then 994 7^ and we obtain 
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fJ-(x) = e (34) 



When combined with Eg. fl25|) , Eqs.(|33]) and alternatively give both an explicit solution to the 
problem and an existence condition, since a solution H (x) - and hence an integrating factor of the 
form /j,(x,y') - exists if the integrand in Eq.(|33|) or Eq.(|3~4|) only depends on x. A 



2.2.3 Proof of Lemma 3 

We start from Eq.(p3|) by considering the expression 

T -<f> G xy {x,y) + G yy {x,y) y' 

y ~ FA^) (35) 

and develop the proof below splitting the problem into different cases. For each case we show how 
to find J-(x,y') satisfying Eq.(|25|), T will then lead to the required integrating factor when, in 

addition to the conditions explained below, the existence conditions for M (x) mentioned in the 
previous subsection are satisfied. 

Case A : G xy /G yy depends on y 

To determine whether this is the case, we cannot just analyze the ratio G xy /G yy itself since it 
is unknown. However, from Eq.(|35|), in this case the factors of T depending on y will also depend 
on y' , and this condition can be formulated as 

When this inequation holds, we determine F y /(x,y') up to a factor depending on x, that is, the 
required F(x, y'), as the reciprocal of the factors of T which depend on y' but not y. A 

Example: Kamke's ODE 226 

This ODE is presented in Kamke's book already in exact form, so we start by rewriting it in explicit 
form as 

y" = x2yV ' t ^ (37) 

y 



x(xy' + 2y) , , 

T = v y ~ — ^ (38) 



We determine T (Eq.(|35|)) as 



The only factor of T containing y is: 



xy' + 2y (39) 
and since this also depends on y' , T(x,y') is given by 
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y') = y' 



(40) 



Case B : either G xy = or G yy = 

When the expression formed by all the factors of T containing y does not contain y' , in Eq,(|36|) 

r) X 

we will have -§^{~y) = 0, and it is impossible to determine a priori whether one of the functions 
{G xy , Gyy} is zero, or alternatively their ratio does not depend on y. We then proceed by assuming 

the former, build an expression for J-(x, y') as in Case A, and check for the existence of M (x) as 

explained in the previous subsection. If A* (x) exists, the problem is solved; otherwise we proceed 
as follows. 

Case C : G xy /G yy ^ and does not depend on y 

In this case, neither G xy nor G yy is zero and their ratio is a function of just x, so that 



Gxy = vi(x) w(x,y) 

G y y = v 2 (x) w(x,y) (41) 
for some unknown functions vi(x) and V2{x). Eq.(|35[) is then given by 



T — w(x, y) <*%+y) (42) 

for some function w(x,y), which is made up of the factors of T depending on y and not on y' . To 
determine F y /(x,y') up to a factor depending on x, we need to determine the ratio v\{x)/v2{x). 
For this purpose, from Eq.(||) we build a PDE for G y (x, y), 

G X y = 7— r Gy y (43) 

y V2(X) yy 



The general solution of Eq. ( |43|) is 

G y = G(y + p(x)) (44) 
where Q is an arbitrary function of its argument and for convenience we introduced 

p'(x) = v\(x)/v2{x) (45) 
We now determine p'(x) as follows. Taking into account Eq. (|4l|) , 

V2(x)w(x,y) = Q\y + p(x)) (46) 
By taking the ratio between this expression and its derivative w.r.t y we obtain 

oy Q'{y + p{x)) 
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that is, a function of y + p(x) only, which we can determine since we know w(x, y). If Tt' 7^ 0, p'(x) 
is given by 



H x w xv w - w x w. 
? 

In summary, the conditions for this case are 



// \ n x w xy w-w x Wy 
P\ x ) = — = ~T~7. 172 ( 48 ) 



TCy WyyW ~ Wy 



T » #0 - w(-t) =0 - ^ [a{vi) * ' ik H "'^ (49) 

and then, from Eq.(f42|), J-(x,y') is given by 

t( >\ (p' + y') w / Kn N 

Fl x , y ) = — y — ^ ' 

where at this point T, w(x,y) and p'{x) are all known. A 
Example: Kamke's ODE 136. 

We begin by writing the ODE in explicit form as 

y" = ^ (51) 

x-y 

This example is interesting since the standard search for point symmetries is made difficult by the 
presence of an arbitrary function of y' . T (Eq.(pq)) is determined as 



(x - y) 2 

and w(x,y) as 



(52) 



Then H{y + p{x)) (Eq.(g7])) becomes 



w(x,y) = - (53) 

(x - y) 2 



n = (54) 

x-y 



p'{x) = -1 (55) 



^ = Wf (56) 

Case D : Ti = 

We now discuss how to obtain p'{x) when TL'{y +p(x)) = 0. We consider first the case in which 
Ti = 0. Then, Q" = and the condition for this case is 



and hence, from Eq.(El), p'(x) is 



so from Eq.(|50|): 
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T y = (57) 

Recalling Eq, (|44]) , G is given by 

G(x, y) =d(y + p{x)f + C 2 (y + p(x)) + g(x) (58) 
for some function g{x) and some constants C\, C 2 . From Eq.(f23"l), $>(x,y,y') takes the form 



M n F x (x,y')+g'(x) + (2d (y + p(x)) + C 2 )(y> + p'{x)) 

*( X >V>V) = F^7) (59) 

We now determine p'(x) as follows. First, from the knowledge of T and <1? we build the two explicit 
expressions: 

1 F i 

A = — = (60) 

T 2d (y'+p'(x)) { ' 



and 



_ <S>{x,y,y') F x +g'(x) C 2 

* = ~^T~ - y= 2 Cl WTpW) +P(X) + 2d (61) 
From Eq,(|60|) and Eq.(61) A and Vl/ are related by: 

§- x Uy' + v'(x)) a) + ^7 (iv' + p'O*)) = p(x) + || (62) 

where the only unknowns are p(x), C±, and C 2 . By differentiating the equation above w.r.t y' and 
x we obtain two equations where the only unknown is p'(x): 



k y ,p"{x) + (A xy , + y y , y ,)(y' +p'(x)) + A X + 2^! y , = (63) 

A p"'{x) + (A x , + Vy' X )(y' + P '{x)) + (A, + %>)p"(x) + V X = p'{x) (64) 

from where we obtain p'(x) by solving a linear algebraic equation built by eliminating p"(x) between 
Eq.(||) and Eq.@)|. Also, as a shortcut, if (A xy > + ^>y>y>)/ A y > depends on y' , then we can build a 
linear algebraic equation for p'(x) by solving for p"(x) in Eq. (|63|) and differentiating w.r.t. y'. A 

Remark 

If Eq.(|63|) depends neither on p'(x) nor on p"(x) this scheme will not succeed. However, in 
that case the original ODE is actually linear and given by Eq. (p8|) . To see this, we set to zero the 
coefficients of p'(x) and p"(x) in Eq.(|63|), obtaining: 



Ay> = k X yl + ^lylyl = A X + 2 y> = (65) 



8 From Eq.( |60|) , A 7^ 0, so that Eq. (pi]) always depends on p"'(x), and solving Eq.(p3|) for p"(x) and substituting 
twice into Eq.(|34D will lead to the desired equation for p'(x). If Eq.(^) depends on p'{x) but not on p"(x), then 
Eq.(p3|) itself is already a linear algebraic equation for p'(x). 
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which implies that A is a function of x only, and then 

Vyy = (66) 

If we now rewrite F{x, y') as 

F(x, y') = Z(x, y') - g(x) - A(y> + p'f C x (67) 



and introduce this expression in Eq.(pO|), we obtain Z y i = 0; similarly, using this result, Eq. ( |6l|) , 
Eq.(|66|) and Eq.(|67|) we obtain Z x = 0. Hence, Z is a constant, and taking into account Eq.(|67|) 
and Eq. (|59|) , the ODE which led us to this case is just a non- homogeneous linear ODE of the form 



(y + p)" + (A'(y + p)' - 2(y + p) - C 2 /d)/2A = (68) 
whose homogeneous part does not depend on p(x): 

2A{x) A{x) 



and as mentioned, it is the same as Eq.(p8[). 
Example: Kamke's ODE 66. 
This ODE is given by 

y" = a(c + bx + y)(y' 2 + lf 12 (70) 
Proceeding as in Case A, we determine T, w(x,y), and H(y + p(x)) as 

T = a(V 2 + l) 3/2 ; w(x,y) = l; H = (71) 

From the last equation we realize that we are in Case D. We determine A and ^ (Eqs. (|60|), (|6l|)) 
as: 



A 



( 3/2 

a 



(^ + i)' 

= c + bx (72) 



We then build Eq.fl62|) for this ODE: 



(/> + !) 

Differentiating w.r.t. y' leads to Eq.(|63|): 

_ 3 P"{x)y' 
" ( y ' 2 + l) 5/Z a 



+c + bx = p{x) + ^- (73) 



3/2 ' ' J 2C 

a 



— ~ " ( 74 ) 
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from which it follows that p"(x) = 0. Using this in Eq.(|64]) we obtain: 

p'(x) = b (75) 

after which Eq,(|50|) becomes 



a (V 2 + l) 

Case E : W = and H + 

In this case TC(y+p(x)) = Q" /Q' = Ci, so Q' is an exponential function of its argument (y+p(x)) 
and hence from Eq. 



G(x, y) = C 2 e^ x ^ + (y + p(x))C 3 + g{x) (77) 
for some constants C 2 , C3 and some function g(x). In this case one of the conditions to be satisfied 



is 



T y = constant 7^ (78) 

and $(x, y, y') will be of the form 

FAx. y') + g'(x) + (c 2 C\e(«+'<*» c > + C\) {y' + p'(x)) 

= Sv(Z7) " (79) 

Taking advantage of the fact that we explicitly know C%, we build a first expression for p' by 
dividing C\e yCl by T: 

Fyl 

A = ~C 2 e^)C, (y' + p '{ x )) (80) 
We obtain a second expression for p' by multiplying $ by A and subtracting C\e ClV 



C 2 eP( x ) Cl \ y' +p'{x) 
As in Case D, A and ^ are related by 

^ Uy' + a) + (y' + A*)) p'{x)Ac t + ^[iv'+ A*)) *J = C2 % Cl ( 82 ) 

where the only unknowns are C 2 , C3 and p(x). Differentiating Eq.(^) with respect to y' we have 

L"(x) + p'(x) 2 C^j A y , + p'{x) (y'ky'Ci + Ad + A xy , + J ^ 

+21y + A;, + y'A^, + j/'*^, = 

The problem now is that, due to the exponential on the RHS of Eq.(|82|), differently from Case D, 
we are not able to obtain a second expression for p'(x) by differentiating w.r.t x. The alternative 
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we have found can be summarized as follows. We first note that if A y i = 0, Eq.(p3|) is already a 
linear algebraic equation^ for p', so that we are only worried with the case Ay ^ 0. With this in 
mind, we divide Eq.(|83|) by K y > and, i/the resulting expression depends on y', we directly obtain a 
linear algebraic equation in p'(x) by just differentiating w.r.t y'. A 

Exampl(^\ 

S( X y + m -2 +e ») 

y x z + y — 1 

We determine T, w(x,y), and TC(y + p(x)) as 

2/(xy' + l)e y . 
T = - r ^- — - — -; w(x,y) = e y ; H = l 85 
y'x 1 + y'-l 

From the last equation we know that we are in Case E. We then determine A and \E' as in Eqs. 
© and (P): 



A _ y'x 2 + y'-l 
y'(xy' + 1) 

* = -2 (86) 

Now, we build Eq.©: 



/ , -i P + P +2/ m- ( & + — J+2sp-2J = — — 87 

xy' + l\\ xy' + lj\ y'J J C 2 e v 



and, differentiating w.r.t. y', we obtain Eq. (|83D : 

./2 



W + l -(x 3 + x)y' 2 / „ >2 n 2?/-l-2x + xy / , , 

Proceeding as explained, dividing by K y i and differentiating w.r.t. y', we have 
9 ( ,2 2 y' - 1 - 2 x + xy' \ , . _ 

' A ; . A r—, rr—, ,2, (V - 1) = (89) 



dy' \ (xy' + 1) (2 xy' + 1 - (x 3 + 



xjy 



We can see this by assuming that A y i = and that Eq. (|83[) does not contain p , and then arriving at a contradiction 
as follows. We first set the coefficients of p in Eq.fcSI) to zero, arriving at 



= Ci A + = 21v + A * + ^y'v'V' ( A ) 

Eliminating ^ y ' y ' gives 

2¥„/ = & A y' - A x 
Differentiating the expression above w.r.t y and since A y i — 0, we have 

2 ^y'y' =Ci A 
Finally, using Eq.(A), = A, contradicting F y > 7^ 0. 

10 There are no examples of this type in all of Kamke's set of non-linear second order ODEs. 
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Solving for p'(x) gives p'(x) = 1/x, from which Eq.(|50|) becomes: 



K*,v) = {v--) + 00) 

Case F 

The final branch occurs when Eq.(^) divided by A„/ does not depend on y' (so that we will not 
be able to differentiate w.r.t y'). In this case we can build a linear algebraic equation for p'(x) as 
follows. Let us introduce the label f3(x,p',p") for Eq. (j33|) divided by A y i, so that Eq. (|83|) becomes: 

A y ,(x,y>) P(x,p',p") = (91) 

Since we obtained Eq.(p3[) by differentiating Eq.(p2|) with respect to y', Eq.(|82[) can be written in 
terms of /3 by integrating Eq.(|9l|) with respect to y'\ 

A(x,y')p(x,p',p") +1 (x,p',p") = c Jl )Cl (92) 

where ^l{x,p' ,p") is the constant of integration, and can be determined explicitly in terms of x, 
p' and p" by comparing Eq.(^) with Eq,(jS2j). Taking into account that (3(x,p',p") = 0, Eq.(|92]) 
reduces to: 

We can remove the unknowns C2 and C3 after multiplying Eq.(|93|) by e p ( x ' Cl , differentiating with 
respect to x, and then dividing once again by e p ^ Cl . We now have our second equation for p', 
which we can build explicitly in terms of p', since we know j(x,p' ,p") and C\\ 

^L + Cl p'j = (94) 

Eliminating the derivatives of p' between Eq,(|9l|) and Eq.(94) leads to a linear algebraic equation 
in p' . Once we have p' , the determination of J-(x, y') follows directly from Eq.([50|). A 



2.3 Integrating factors of the form fi(y, y') 

From Eq,(|ll|), the ODE family admitting an integrating factor of the form n(y,y') is given by 

*" = -£( G » + |/"*')-T (95) 

where fj,(y, y') and G(x, y) are arbitrary functions of their arguments. For this ODE family, it 
would be possible to develop an analysis and split the problem into cases as done in the previous 
section for the case n(x,y'). However, it is straightforward to notice that under the transformation 



y{x) — > x, x — > y(x), Eq.(^) transforms into an ODE of the form Eq.(22) with integrating factor 
fx(x, y'~ 1 )/y' 2 . It follows that an integrating factor for any member of the ODE family above can be 
found by merely changing variables in the given ODE and calculating the corresponding integrating 
factor of the form fj,(x,y'). 
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Example: 



/2 

y 

y" h sin(x) y 1 y + cos(x) y 2 = (96) 



Changing variables y(x) — ► x, x — > y(x) we obtain 



y" H sin(y) y x — cos{y)x y' = (97) 

Using the algorithm outlined in the previous section, an integrating factor of the form fi(x, y') for 
Eq.(p7[) is given by 

4- (98) 



y' x 



from where an integrating factor of the form fi(y,y') for Eq.(|96|) is 1/y, leading to the first integral 



sin(x) y + - + d = 0, (99) 



which is a first order ODE of Bernoulli type. The solution to Eq.(|96|) then follows directly. This 
example is interesting since from [|j Eq.(96) has no point symmetries. 



3 Integrating factors and symmetries 

Besides the formulas for integrating factors of the form /i(x,y), the main result presented in this 
paper is a systematic algorithm for the determination of integrating factors of the form /j,(x, y') and 
/i(y,y') without solving any auxiliary differential equations or performing differential Groebner basis 
calculations, and these last two facts constitute the relevant point. Nonetheless, it is interesting 
to briefly compare the standard integrating factor (fi) and symmetry approaches, so as to have an 
insight of how complementary these methods can be in practice. 

To start with, both methods tackle an n th order ODE by looking for solutions to a linear 

n th 

order determining PDE in n + 1 variables. Any given ODE has infinitely many integrating factors 
and symmetries. When many solutions to these determining PDEs are found, both approaches can, 
in principle, give a multiple reduction of order. 

In the case of integrating factors there is one unknown function, while for symmetries there is a 
pair of infinitesimals to be found. But symmetries are defined up to an arbitrary function, so that 
we can always take one of these infinitesimals equal to zero0; hence we are facing approaches of 
equivalent levels of difficulty and actually of equivalent solving power. 

Also valid for both approaches is the fact that, unless some restrictions are introduced on the 
functional dependence of (i or the infinitesimals, there is no hope that the corresponding determining 
PDEs will be easier to solve than the original ODE. In the case of symmetries, it is usual to restrict 
the problem to ODEs having point symmetries, that is, to consider infinitesimals depending only 

11 Symmetries [£,(x, y, ..y^ n ~ 1 ' > ),r)(x, y, -y^ n ~ 1 ')] of an n th order ODE can always be rewritten as [G, (G — + rj\, 
where G(x, y, ..j/" -1 ') is an arbitrary function (for first order ODEs, y' must replaced by the right-hand-side of the 
ODE). Choosing G — the symmetry acquires the form [0, fj\ 
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on x and y. The restriction to the integrating factors here discussed is similar: we considered /i's 
depending on only two variables. 

At this point it can be seen that the two approaches are complementary: the determining 
PDEs for {i and for the symmetries are differential, so that even using identical restrictions on the 
functional dependence of fx and the infinitesimals, problems which may be untr actable using one 
approach may be easy or even trivial using the other one. 

As an example of this, consider Kamke's ODE 6.37 

y" + 2 y y' + f(x) (y' + y 2 ) - g(x) = (100) 

For arbitrary f(x) and g(x), this ODE has an integrating factor depending only on x, easily de- 
termined using the algorithms presented. Now, for non-constant f(x) and g(x), this ODE has no 
point symmetries, that is, no infinitesimals of the form [£(x,y), n(x, y)}, except for the particular 
case in which g(x) can be expressed in terms of f(x) as inP| 

jn 3 jy/ C 2 exp^-3/2 J f(x)dx 

4 (2d + J exp (-1/2 J f(x)dx\ dx 

Furthermore, this ODE does not have non-trivial symmetries of the form [£(x,y'), rj(x,y')] either, 
and for symmetries of the form [£(y, y'), r/(y, y')] the determining PDE does not split into a system. 
Another ODE example of this type is found in a paper by || (1988): 

y" - — - g{x) P y p y' - g'y p+1 = (102) 
y 

In that work it is shown that for constant p, the ODE above only has point symmetries for very 
restricted forms of g(x). For instance, Eq.(|96|) is a particular case of the ODE above and has no 



g (x) = J — + -^J- + t ^ 1 L ? (101) 
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point symmetries. On the other hand, for arbitrary g(x), Eq.(102) has an obvious integrating factor 



depending on only one variable: 1/y, leading to a first integral of Bernoulli type: 

V -- g (x)yV + C l = U (103) 

so that the whole family Eq. (|102| ) is integrable by quadratures. 

We note that Eq. (|l00| ) and Eq. (|l02| ) are respectively particular cases of the general reducible 
ODEs having integrating factors of the form /u(x): 

y" = -i^±^-y' - (104) 

where [i{x) and G(x,y) are arbitrary; and n(y): 
12 We are considering here ODEs of order greater than one. 

13 To determine g(x) in terms of f(x) we used the standard form Maple package by Reid and Wittkopf complemented 
with some basic calculations. 



17 



My) Kv) 



In turn, these are very simple cases if compared with the general ODE families Eq.(|23|) and Eq.(|95|), 
respectively having integrating factors of the forms fj,(x,y') and n(y,y'), and which can be system- 
atically reduced in order using the algorithms here presented. 

It is then natural to conclude that the integrating factor and the symmetry approaches are 
useful for solving different types of ODEs, and can be viewed as equivalently powerful and general, 
and in practice complementary. Moreover, if for a given ODE, an integrating factor and a symmetry 
are known, in principle one can combine this information to build two first integrals and reduce 
the order by two at once (see for instance 



4 Tests 

After plugging the reducible-ODE scheme here presented into the ODEtools package |i~0f| , we tested 
the scheme and routines using Kamke's non-linear 246 second order ODE examples^ The purpose 
was to confirm the correctness of the returned results and to determine which of these ODEs have 
integrating factors of the form /i(x,y), ji{x,y') or fj,(y,y'). The test consisted of determining fx and 
testing the exactness condition Eq.(||). 

In addition, we ran a comparison of performances in solving a subset of Kamke's examples 
having integrating factors of the forms fi(x,y') or /z(y,y'), using different computer algebra ODE- 
solvers (Maple, Mathematica, MuPAD and the Reduce package Convode). The idea was to situate 
the new scheme in the framework of a sample of relevant packages presently available. 

To run the comparison of performances, the first step was to classify Kamke's ODEs into: 
missing x, missing y, exact and reducible, where the latter refers to ODEs having integrating factors 
of the forms fj,(x,y') or fi(y,y'). ODEs missing variables were not included in the test since they 
can be seen as first order ODEs in disguised form, and as such they are not the main target of the 
algorithm being presented. The classification we obtained for these 246 ODEs is as follows 



Classification 


ODE numbers as in Kamke's book 


99 ODEs are missing x 
or missing y 


1, 2, 4, 7, 10, 12, 14, 17, 21, 22, 23, 24, 25, 26, 28, 30, 31, 32, 40, 42, 43, 45, 46, 47, 48, 
49, 50, 54, 56, 60, 61, 62, 63, 64, 65, 67, 71, 72, 81, 89, 104, 107, 109, 110, 111, 113, 
117, 118, 119, 120, 124, 125, 126, 127, 128, 130, 132, 137, 138, 140, 141, 143, 146, 150, 
151, 153, 154, 155, 157, 158, 159, 160, 162, 163, 164, 165, 168, 188, 191, 192, 197, 200, 
201, 202, 209, 210, 213, 214, 218, 220, 222, 223, 224, 232, 234, 236, 237, 243, 246 


13 are in exact form 


36, 42, 78, 107, 108, 109, 133, 169, 170, 178, 226, 231, 235 


40 ODEs are reducible 
with integrating factor 
fi(x,y') or fi(y,y') and 
missing x or y 


1, 2, 4, 7, 10, 12, 14, 17, 40, 42, 50, 56, 64, 65, 81, 89, 104, 107, 109, 110, 111, 125, 126, 
137, 138, 150, 154, 155, 157, 164, 168, 188, 191, 192, 209, 210, 214, 218, 220, 222, 236 


28 ODEs are reducible 
and not missing x or y 


36, 37, 51, 66, 78, 97, 108, 123, 133, 134, 135, 136, 166, 169, 173, 174, 175, 176, 178, 
179, 193, 196, 203, 204, 206, 215, 226, 235 



Table 1. Missing variables, exact and reducible Kamke's 246 second order non- linear ODEs. 



Kamke's ODEs 6.247 to 6.249 cannot be made explicit and are then excluded from the tests. 
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For our purposes, the interesting subset is the one comprised of the 28 ODEs not already missing 
variables. The results we obtained using the aforementioned computer algebra ODE-solvers^] are 
summarized as follows^: 



Kamke's ODE numbers 




Convode 


Mathematica 3.0 


MuPAD 1.3 


ODEtools 


Solved: 


51, 166, 173, 
174, 175, 176, 
179. 


78, 97, 108, 166, 
169, 173, 174, 175, 
176, 178, 179, 206. 


78, 97, 108, 133, 
166, 169, 173, 174, 
175, 176, 179, 


51, 78, 97, 108, 133, 134, 135, 136, 
166, 169, 173, 174, 175, 176, 178, 
179, 193, 196, 203, 204, 206, 215. 


Totals: 


7 


12 


11 


22 


Reduced: 








36, 37, 66, 123, 226, 235. 


Totals: 
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Table 2. Performances in solving 28 Kamke's ODEs having an integrating factor fi(x, y) or [i(y, y ) 



As shown above, while the scheme here presented is finding first integrals in all the 28 ODE 
examples, opening the way to solve 22 of them to the end, the next scores are only 12 and 11 ODEs, 
respectively solved by Mathematica 3.0 and MuPAD 1.3. 

Concerning the six reductions of order returned by odsolve, it must be said that neither MuPAD 
nor Mathematica provide a way to convey them, so that perhaps their ODE-solvers are obtaining 
first integrals for these cases but the routines are giving up when they cannot solve the problem to 
the end. 

5 Conclusions 

In connection with second order ODEs, this paper presented a systematic method for determining 
the existence of integrating factors and their explicit form, when they have the forms fj,(x, y), [i(x, y') 
and fi(y, y'). The scheme is new, as far as we know, and its implementation in the framework of the 
computer algebra package ODEtools has proven to be a valuable tool. Actually, the implementation 
of the scheme solves ODEs not solved by using standard or symmetry methods (see sec. |||) or some 
other relevant and popular computer algebra ODE-solvers (see sec. ||). 

Furthermore, the algorithms presented involve only very simple operations and do not require 
solving auxiliary differential equations, except in one branch of the fi(x, y) problem. So, even for 
examples where other methods also work, for instance by solving the related PDE system Eqs.(^) 
and (|7|) using ansatzes and differential Groebner basis techniques, the method here presented can 
return answers faster and avoiding potential explosions of memoryp]. 

On the other hand, we have restricted the problem to the universe of second order ODEs having 
integrating factors depending only on two variables while packages as CONLAW (in REDUCE) can 
try and in some cases solve the PDE system Eqs.(|6|) and (|7|) by using more varied ansatzes for fi. 

15 Maple R4 is not present in the table since it is not solving any of these 28 ODEs. This situation is being resolved 
in the upcoming Maple R5, where the ODEtools routines are included in the Maple library, and the previous ODE- 
solver was replaced by odsolve. However, the scheme here presented was not ready when the development library 
was closed; the reducible scheme implemented in Maple R5 is able to determine, when they exist, integrating factors 
only of the form fj,(y'). 

16 Some of these 28 ODEs are given in Kamke in exact form and hence they can be easily reduced after performing 
a check for exactness; before running the tests all these ODEs were rewritten in explicit form by isolating y". 

1 Explosions of memory may happen when calculating all the integrability conditions involved at each step in the 
differential Groebner basis approach. 
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A natural extension of this work would be to develop a scheme for building integrating factors of 
restricted but more general forms, now for higher order ODEs. We are presently working on these 
possible extensions^, and expect to succeed in obtaining reportable results in the near future. 
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Appendix A 



We display here both the integrating factors obtained for the 28 Kamke's ODEs used in the tests 
(see sec. ffl) and the "case" corresponding to each ODE when using just the algorithms for fi(x, y') 



or n(y, As explained in sec. 2.2.5 , the algorithm presented is subdivided into different cases: 

A, B, C, D, E and F, and case B is always either A or C. 



Integrating factor 



Kamke's book ODE-number 



Case 



o J' f(x)dx 



36 



D 



37 



A 



b+y' 



51, 166, 169, 173, 175, 176, 179, 196, 203, 204, 206, 
215 



C 



(i+y' 2 Y 



66 



D 



78 



D 



97 



A 



108 



D 



123 



A 



133 



C 



(i+v')v' 



v'-i 
(i+j/')(nV 2 ) 



134 



c 



135 



C 



y'-i 
h (v') 



136 



C 



2 xy' — l 



(i + y)' 1 



174 



C 



178 



C 



y 1 (1 + 2 yy') 



193 



C 



y 



226 



A 



235 



C 



Integrating factors for Kamke's ODEs which are reducible and not missing x or y. 



20 We note that for non-linear ODEs these two algorithms work as well when /i y < = 0, but in practice these very 
simple examples are covered by the algorithm for /i(x,y) presented in sec. 2.1. 
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